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Abstract 


A search for pair production of neutral color-octet weak-triplet scalar particles (0°) is 
performed in processes where one 0° decays to a pair of b quark jets and the other to 
a Z boson plus a jet, with the Z boson decaying to a pair of electrons or muons. The 
search is performed with data collected by the CMS experiment at the CERN LHC 
corresponding to an integrated luminosity of 19.7 fb -1 of proton-proton collisions at 
\/s = 8 TeV. The number of observed events is found to be in agreement with the 
standard model predictions. The 95% confidence level upper limit on the product of 
the cross section and branching fraction is obtained as a function of the 0° mass. The 
95% confidence level lower bounds on the 0° mass are found to be 623 and 426 GeV, 
for two different octo-triplet theoretical scenarios. These are the first direct experi¬ 
mental bounds on particles predicted by the octo-triplet model. 
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1 Introduction 

The CERN LHC is well suited to searches for new colored particles. These searches have tar¬ 
geted particles including vector-like quarks (color triplet spin-1/2 particles), colorons (color 
octet spin-1 particles), gluinos (color octet Majorana fermions), and color-octet weak-singlet 
scalars [1 -6|. With the recent discovery of a Higgs boson [ZHU, so far consistent with being 
a neutral scalar colorless particle [l lOl [ill , it is also interesting to look for additional scalars, 
in this case with color charge, with masses close to the electroweak scale. Particles with this 
combination of spin, charge, and color quantum numbers have not been thoroughly sought at 
collider experiments. Neutral color-octet spin-0 particles (0°), for example, emerge in the octo- 
triplet model ltT2l . which includes three particles (0 1 ,0°,0 ) that transform as a color octet 
and weak triplet under SU( 3) c x SU(2) w . The 0° may arise from the sector responsible for 
breaking extended gauge symmetry associated with colorons 10314151 . Octo-triplets may also 
be fermion-antifermion bound states HT6H T7L or elementary particles from non-minimal grand 
unified theory. To date, no direct experimental bound has been set for particles predicted by 
the octo-triplet model. 

Octo-triplet particles would be produced in pairs at the LHC either through quark-antiquark 
annihilation or gluon-gluon interactions, with the former mediated by coloron (G') particles or 
by gluons (g). If a G' is produced, it can decay to 0°©°, ©°(pi, or two quarks, where (pi is a 
color singlet from a general renormalizable coloron model Ifl5l . In the present analysis we are 
searching for 0° produced in pairs. The (pi mass (m^) is a free parameter and has an indirect 
impact on the 0° pair production cross section, via the G' width. 

Unlike color-octet weak-singlet particles, octo-triplet particles do not decay to gluons. How¬ 
ever, with a 0° in the loop, one-loop decays to an electroweak boson and a gluon are allowed. 
There is a potentially similar rate of decays to standard model (SM) quark pairs via mixing 
through a vector-like quark, which can be as heavy as a fewTeV. These heavy vector-like 
quarks appear in many extensions of the SM, such as composite Higgs models fl8l and little 
Higgs models Il9l . In the former, the Higgs boson may be a bound state of a top quark and the 
heavy vector-like quark, where the binding is provided by the spin-1 G' coloron. 

In this paper we assume that one 0° decays to two b quark jets and the other to a Z boson plus 
a jet originating from a gluon, with Z —> ££, where £ is an electron or muon. With no existing 
constraints on the ratio of 0° decays to quarks versus those decaying to Zg, we have studied a 
final state that includes both decay modes. Requiring Z -^ ££ provides a strong experimental 
signature for event selection and background suppression. The leading order (LO) Feynman 
diagrams for 0° pair production by quark-antiquark annihilation and gluon-gluon interaction, 
and the decays explored in this analysis, are shown in Fig. |T| 
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Figure 1: Leading order Feynman diagrams for 0° pair production by quark-anti quark anni¬ 
hilation, through an intermediate G' or gluon (left), and the gluon-gluon interaction (right). In 
addition to the s-channel gluon process shown in the figure, the gluon-gluon interaction can 
also proceed directly through gg —> 0°0° or a f-channel 0° exchange. The decays of the 0° 
and Z boson are also included, where the 0° decays are described in the main text. 
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1 Introduction 


The signal production cross section for the octo-triplet model fT2j is evaluated at LO with 
the MadGraph v4.5.1 program l20l . Figure [2] shows the values of the cross section in the 
2-dimensional plane of G' mass (me,') and 0° mass (m 0 o). The masses of 0°, G', and (pi are free 
parameters in general, but we use = 125 GeV + m & o/3 as a benchmark. We explore two 
distinct mass scenarios: niQi = 2.3m 0 o and niQi = 5m 0 o. The first example corresponds to a 
scenario in which the signal cross section is enhanced because of the additional contributions 
from colorons (see Fig. [2j|, while the second relation corresponds to a scenario in which the 
contribution from colorons is negligible. For both mass relations we use 0° masses, ranging 
from 200 to 900 GeV, which are large enough to allow 0° decay to on-shell Z bosons, but small 
enough to be within the sensitivity of this search. The first mass relation, me = 2.3»z@o, results 
in 0° pair production cross sections from 125 pb down to 2.23x10 2 pb, for the range of 0° 
masses considered. This mass relation also sets the G' mass sufficiently above 0° pair produc¬ 
tion threshold, but small enough compared to m 0 o that production through G' is considerable, 
at 49 to 86% of the total, depending on m 0 o. In this case, quark-antiquark annihilation dom¬ 
inates, while gluon-gluon interaction production decreases with increasing 0° (and G') mass, 
because of the differences in the parton distribution functions (PDF) of quarks and gluons 112X1 . 
The second mass relation considered, m 0 / = 5m@ o, results in 0° pair production cross sections 
from 63.3 pb down to 2.62 x 10 3 pb, for the range of 0° masses considered, which are a factor 
of two smaller than for the first mass relation. This corresponds to a region where gluon-gluon 
interactions dominate and production through G' is just 0.6 to 3.6% of the total, depending on 
the 0° mass. The cross sections discussed above agree with those calculated in Ref. 1221 where 
a similar model is considered. 
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Figure 2: Cross sections for 0° pair production, as a function of the coloron G' and the octo- 
triplet 0° masses. The dotted lines show a few examples of assumed linear relationship be¬ 
tween the two masses. 

The decay modes 0° —» Zg and 0° —» yg are invariant under SU (2)w- However, in the La- 
grangian, the coefficients of these modes turn out to be much smaller than one. This feature mo¬ 
tivates consideration of other decay modes, which can be represented by higher-dimensional 
operators, and therefore are usually considered as negligible with respect to dominant decay 
modes. In particular, qcj decay modes could be as large as the electroweak boson plus gluon 
modes Ifl2l . In the 0° mass range we consider, the dominant fermion decay mode is either 
bb or tt depending on the values of coefficients appearing in different terms of the 0° -X qq 
Lagrangian. As a benchmark, we assign half of the branching fraction to 0° —> bb, leaving 
the other half to the 0° —> Zg and 0° —> 7 g decay modes (i.e. B(©° —> tt) = 0). The ratio of 
branching fractions £>(©° -X Zg )/B(©° —> 7 g) is set by the mass of the 0°. For the range of 
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0° masses explored, £>(©° —> Zg) is to 38% and 6(0° —> 7 g) accounts for the remaining to 
17%. The total cross section times branching fraction for the final state studied in this analysis 
ranges from 2.77 pb down to 5.74 x 10 4 pb for the mass relation tuqj = 2.3m qo, and 1.40 pb 
down to 6.73 x 1CV 5 pb for the mass relation Mqi = 5m@o. 


2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diame¬ 
ter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip 
tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron 
calorimeter, each composed of a barrel and two endcap sections. Muons are measured in gas- 
ionization detectors embedded in the steel flux-return yoke outside the solenoid. Extensive 
forward calorimetry complements the coverage provided by the barrel and endcap detectors. 
A more detailed description of the CMS detector, together with a definition of the coordinate 
system used and the relevant kinematic variables, can be found in Ref. Il23l . 


3 Event selection 

This analysis uses proton-proton collision data at a center-of-mass energy of 8 TeV collected 
with dilepton (ee, pp, and ep) triggers with transverse momentum pj thresholds of 17 and 
8 GeV for the two lepton candidates. The ep data set is used to estimate tt background. The 
data sample corresponds to an integrated luminosity of 19.7 fb 1 . We select events with a final 
state containing a Z boson and a jet (Z + jet), together with two b quark jets (b jet pair), where 
the Z boson decays to a pair of electrons or muons. 

Electrons are reconstructed using selection criteria that take into account radiated photons 1124] . 
Candidates are required to have pj > 20 GeV and \rj\ < 2.5 where // is the pseudorapidity. 
Electrons from photon conversions are rejected. An isolation condition is imposed within a 
cone of A R = V(A ?/) 2 + (A cp ) 2 < 0.3 around the electron. The sum of the pj of charged 
particles and transverse energy Ej = E sin 9 of neutral particles within this cone is corrected 
in each event for energy deposits due to additional interactions within beam bunch crossings 
(pileup). The corrected sum must be less than 15% of the electron pj. 

Muons are reconstructed using selection criteria based on quantities measured in the tracker 
and muon sub-detectors. Candidates are required to have pj > 20 GeV and |?/| < 2.4. The 
track associated with the muon candidate is required to have hits in the pixel, silicon strip, and 
muon systems. An isolation condition requires that the sum of the pj of charged particles and 
Ej of neutral particles within a cone of A R < 0.4 around the muon is less than 12% of the muon 
pj. The pileup correction for muon isolation is similar to the one applied for electrons. 

A particle-flow technique H25l |26 | is used to identify jet constituents, which are input to the 
FASTjET algorithm Il27l [281 for clustering using the anti -kj algorithm H29l with a distance pa¬ 
rameter of 0.5. The jet energy scale (JES) is measured in data with Z/ 7 +jet and dijet events 
130, [3TI . and a correction is applied to both data and simulated samples. The corrected jets 
must have pj > 40 GeV and \p\ < 2.4. Jets within A R = 0.3 from an isolated electron or muon 
as defined above are not counted as part of the 3 jet requirement. Additional requirements, 
based on the energy balance between charged and neutral hadronic energy in the jet, are used 
to reduce contamination from misidentified jets. A jet is considered b tagged, consistent with 
originating from the hadronization of a bottom quark, if it satisfies the "loose operating point" 
requirements of the combined secondary vertex tagger 11521 . The tagging efficiency is about 
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70-90% and mistagging rate is about 10-20%, depending on pj and rj of the jets. 

Candidate events must satisfy the following criteria: at least one reconstructed primary vertex 
satisfying |z| < 24 cm and impact parameter less than 2 cm; an opposite-sign, same-flavor 
lepton (electron or muon) pair, sharing a primary vertex and with an invariant mass between 
80 and 100 GeV to form a Z boson candidate; the two leptons separated from jets by A R > 0.5; 
at least two b-tagged jets forming a b jet pair system (sum of the b jet four-momenta); and at 
least one additional jet. If more than two b-tagged jets are present in the event, the jets with the 
largest pj values are selected for the b jet pair system. 

The mass of one 0° candidate is reconstructed from the b jet pair system, and the mass of 
the other 0° candidate is reconstructed from the combination of the Z boson candidate and 
the highest fj jet that is not a part of the b jet pair system. With this prescription, the correct 
combination of the Z boson and the jet in simulated signal samples is chosen about 65 to 80% 
of the time, depending on 7n & o. 

For each hypothetical mass point, we consider signal events in a rectangle that is formed in the 
two-dimensional mass plane defined by the two reconstructed 0° masses. We determine the 
center, length, and width of the rectangle for each mass point by fitting the b jet pair and Z + jet 
mass distributions in the simulated signal samples with a modified Gaussian function that has 
an additional parameter to account for low-mass non-Gaussian tails. The length and width 
of each rectangle are defined as ±3<j sd where cr S(d is the standard deviation from the modified 
Gaussian fit. For example, the rectangle is 240 GeV long and 130 GeV wide for the signal with 
niQi = 1100 GeV and m©o = 478 GeV (Fig. [3] (left)). We choose a 3cr sd wide window to keep more 
signal events in middle-high mass regions, since most of the SM background events appear in 
low-mass region (Fig. [3] (right)). The signal search regions are defined for the electron and 
muon channels separately, although these turn out to be very similar. 



b jet pair mass [GeV] b jet pair mass [GeV] 


Figure 3: Distributions of Z + jet mass versus b jet pair mass for signal events with m G / = 
1100 GeV and m©o = 478 GeV (left) and background events (right), in the electron channel 
(distributions are similar in the muon channel). The open red rectangle on left indicates the 
signal region, as described in the text. 


4 Backgrounds 

The dominant background processes are Z + jets and tt + jets. The contributions from other 
processes, including diboson+jets and multijet, are small. The samples used for all back¬ 
ground processes are generated at tree level with the MadGraph program, interfaced with 
PYTHIA v6.4 |33| for showering and hadronization. The Z + jets background is normalized to 
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the next-to-next-to-leading-order cross section using FEWZ v2.1 Il34l . The tt + jets and tW+jets 
(single top quark) backgrounds are normalized to the next-to-next-to-leading-logarithm cross 
sections |35l |36j|- The diboson+jets backgrounds are normalized to the next-to-leading-order 
(NLO) cross section from the MCFM v5.8 ||37| calculation. Full simulation of the CMS detector 
is implemented using the GEANT4 package [[38] . 

4.1 Z + jets 

More than half of the background events come from Drell-Yan (DY) production in association 
with jets. There are two points in estimating the background that require special attention. 
First, care must be taken in estimating the yield of Z + >3 jets events as the simulation may not 
correctly model the kinematic properties of multijets. Second, the event yield and shape in the 
Z + b-jets process, where the kinematic properties of b jet pairs and the fraction of heavy-flavor 
jets in inclusive Z + jets might also suffer from mis-modeling in the simulation. The following 
paragraphs describe how the data are used to improve on the estimates from simulations. 

To address the first point, the event yield from the simulation is multiplied by a correction factor 
to normalize it to data in a control region. This control region is defined such that the dilepton 
mass is between 80 and 100 GeV and there are at least three jets, none of which is b tagged. The 
correction factor is found to be 0.98 ± 0.12 (0.91 ± 0.12) for the electron (muon) channel. The 
correction is applied to all signal rectangle regions. The 12% uncertainty comes from b tagging, 
JES, and pileup systematic uncertainties in the control region, added in quadrature, while the 
statistical uncertainty is negligible. The Z + jet mass in the Z + >3 jets (one b tag) control region 
is plotted in Fig.[4]after the correction factor has been applied. Agreement between the data and 
simulated sample is observed in the Z + >3 jets (no b tag) and Z + >3 jets (one b tag) control 
regions. 

The second point requires a different approach since control regions that include two or more 
b jets may suffer from signal contamination. We take a two-step approach. First, the simulated 
Z + b-jets events are weighted by the ratio of the k-factor (NLO cross section divided by LO 
cross section) for Z + 2 b-jets to that of Z+2 jets, using MCFM [[391 . The ratios vary from 1.09 to 
2.56 in the b jet pair mass range of 20 GeV to 1.8 TeV. In the simulation, about 20% of Z + jets 
events in the signal region have at least two jets originating from b quarks. In the second 
step, the remaining difference between data and simulation is evaluated as a function of b jet 
pair mass in the non-signal, off-diagonal regions (sidebands) in the b jet pair mass and Z + jet 
mass plane. The uncertainty in the Z+heavy-flavor jets processes is taken from this difference 
or from the uncertainty in the CMS cross section measurement (40] , whichever is larger. The 
uncertainty varies from 20 to 50%. 

The uncertainty in the shape of the other variable that is used to define the signal search regions, 
the Z + jet mass distribution, is estimated by comparing the distribution in the simulated Z + 
jets sample with those of several MadGraph samples that are produced with factorization or 
renormalization scales and matrix element parton shower matching thresholds varied up and 
down by a factor of two. The maximum difference between the nominal distribution and the 
varied distributions is taken as an uncertainty The uncertainty varies from 2 to 55%. 

4.2 tt + jets and tW + jets 

After accounting for the Z + jets background, most of the remaining background events come 
from tt + jets, with a smaller contribution from tW+jets production. 

Both tt + jets and tW+jets processes can yield final states containing an opposite-sign e» pair. 
In this analysis, we are considering only opposite-sign, same-flavor lepton pairs, and thus the 
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Figure 4: The Z + jet mass distribution in the Z + >3 jet (one b tag) control region after the 
appropriate correction factor has been applied, for the electron channel (left) and muon channel 
(right). The panels at the bottom show the ratio of data to background simulation, with the 
band representing the systematic uncertainty including the normalization uncertainty obtained 
from the Z + >3 jet (no b tag) control region. 


data containing an opposite-sign e// pair can be used to calibrate the tt + jets and tW+jets back¬ 
grounds. The overall event yield in the simulation is multiplied by a factor of 1.07 ± 0.13 to 
normalize it to data in the control region containing an e}i pair with the mass between 60 and 
120 GeV (Z efl ) and at least three jets, at least two of which are b tagged. The uncertainty comes 
from b tagging, JES, and pileup systematic uncertainties in the control region, added in quadra¬ 
ture, while the statistical uncertainty is negligible. The b jet pair and Z e/( +jet mass distributions 
in this control region are plotted in Fig. [5] after the normalization factor is applied. 



Figure 5: The Z e „+jet (left) and b jet pair mass (right) distributions in the tt control region after 
the appropriate correction factor is applied. The panels at the bottom show the ratio of data to 
background simulation, with the band representing the systematic uncertainty. 


4.3 Diboson+jets 

A few percent of the background events come from diboson+jets processes. The simulated 
diboson cross sections are found to be in agreement with those measured by CMS ET[ [42 j. 
Therefore we make no correction to the event yields. 

4.4 Multijet events 

The estimation of the multijet background with simulated samples is difficult given that the 
mis-identification rates of jets as electrons or muons are small [2 4 1143 | . Unlike Z -+ it and 
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W W + —> Itvv processes, there are about equal amount of like-sign and opposite-sign lepton 
pairs in the background because the leptons are either misidentified or are non-prompt. In this 
analysis, we consider only opposite-sign, same-flavor lepton pairs, and thus the data contain¬ 
ing like-sign, same-flavor lepton pairs can be used as a control sample to determine the size 
of the contribution from multijet events. After the Z boson mass window selection and the 
requirement of at least three jets, with at least two of them being b tagged, and after subtract¬ 
ing events from all other simulated background processes, 27 (13) dielectron (dimuon) events 
are left in the like-sign dilepton channel. This corresponds to 1.6% (0.7%) of the number of 
opposite-sign dilepton events. These events are distributed mostly in the lower mass region, 
as is the case for other backgrounds. Given the small number of like-sign events in the control 
sample and their mass distribution, the multijet background is found to be negligible. 


5 Systematic uncertainties 

We organize the sources of systematic uncertainties in three categories: sources affecting both 
the signal and background; sources affecting only the background; and sources affecting only 
the signal. A summary of the systematic uncertainties is given in Table [l| 

Sources affecting both the signal and background: 

• b tagging 

The data-to-simulation correction factor for b tagging is varied by its uncertainties. 

The resulting uncertainties in the number of predicted signal and background events 
range from 13 to 25%. 

• Jet energy scale 

The jet energies are varied by the uncertainty in the applied JES correction. The 
resulting event yield uncertainties range from 0.2 to 2.6% for the signal and 4 to 8% 
for the background. 

• Lepton trigger and identification (ID) 

The lepton trigger and ID (including isolation) uncertainty is estimated by vary¬ 
ing the data-to-simulation correction factor by its uncertainty, which is measured in 
dilepton data where at least one of the leptons passes stringent identification criteria. 

The yield uncertainty is about 1%, both for the signal and the background. 

• Pilenp modeling 

An uncertainty in the pileup modeling in the signal and background samples is 
estimated by varying by 5% the total inelastic cross section, as measured at the 
LHC Il44l [45j. The event yield uncertainty is about 1%, somewhat higher for the 
signal than for the background. 

Sources affecting only the background: 

• Z + jet mass shape for Z + jets background 

The uncertainty is evaluated by comparing the varied simulated samples as men¬ 
tioned in Section[4j The shape of the Z + jet mass distribution is much less certain in 
the high-mass tail than in the low-mass region because of the limited number of sim¬ 
ulated events with high masses. The uncertainty in the estimated number of Z + jets 
events varies from 2 to 55%. 

• b jet pair mass shape for Z + jets background 

The uncertainty from this source is evaluated in the off-diagonal sidebands in the 
b jet pair mass and Z+jet mass plane, as mentioned in Section |4j and takes into ac- 
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5 Systematic uncertainties 


count the uncertainty in the CMS cross section measurement of Z + heavy-flavor jets 
processes. The resultant uncertainty in the Z + jets event yield is in the range 20 to 
50% 

• Normalization of Z +jets and tt + jets backgrounds 

The uncertainty in the event yield is 12% for Z + jets and 13% for tt + jets processes 
as mentioned in Section [4] 

• Diboson cross section 

The uncertainty is taken from the underlying CMS diboson cross section measure¬ 
ments fiTl |42| , and implies a 10% uncertainty in the estimated number of diboson 
events. 

Sources affecting only the signal: 

• Initial state radiation (ISR) modeling 

For the signal samples, the uncertainty from this source is estimated from the pj 
distribution of the ©°0° system. The distribution generated using MADGRAPH is 
found to be in agreement with the data for heavy object systems with pj < 120 GeV 
Il46l . Above 120 GeV, an uncertainty is assigned that varies from 5 to 20%, depending 
on the pj. This leads to a 1-5% uncertainty in the yield of signal events for the mass 
ranges that are considered. 

• Integrated luminosity 

The CMS experiment collected data equivalent to 19.7 fb 1 with a 2.6% uncertainty 
Ii47l . This uncertainty applies only to yield of signal events. The background events 
are either normalized in a control region or the normalization is taken from data. 

The systematic uncertainties in event yield due to the lepton energy scale, and the lepton and 
jet energy resolution, are studied and found to be negligible (<0.1%) compared to the other 
sources of uncertainty. 

Table 1: Impact of systematic uncertainties on individual event yields. Ranges show the vari¬ 
ation over the search regions that are considered. Dashes indicate cases where a systematic 
uncertainty is not applied. Sources appearing in more than one process are treated as corre¬ 
lated in the limit setting. 


Source 

Signal [%] 

Z + jets [%] 

tt and tW [%] 

Diboson [%] 

Signal and background 

b tagging 

13-25 

15-16 

13-15 

16-17 

Jet energy scale 

0.2-2.6 

6-8 

4-6 

5-7 

Lepton ID, isolation, trigger 

0.9-1.2 

0.9-1.2 

0.9-1.3 

0.9-1.2 

Pileup modeling 

0.1-1.5 

0.3-0.7 

0.2-1.0 

0.3-0.7 

Background only 

Z + jet mass shape 

— 

2-55 

— 

— 

b jet pair mass shape 

— 

20-50 

— 

— 

Normalization 

— 

12 

13 

— 

Diboson cross section 

— 

— 

— 

7-10 

Signal only 

ISR 

1-5 

— 

— 

— 

Integrated luminosity 

2.6 

— 

— 

— 


For the signal samples, the impact on the event yields due to PDF uncertainties are estimated 
by following the PDF4LHC recommendation [481:49:]. These uncertainties are not used in the 
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limit setting, but are instead included as bands on the theoretical predictions. Since the Z + 
jets and tt + jets backgrounds are normalized in control regions, PDF uncertainties are not 
applied. Similarly the diboson+jets background normalization is taken from the data [41,42j. 
The uncertainty in the signal yield is 7-45%. 


6 Results 


The numbers of observed data events in each signal region after all of the selection require¬ 
ments are applied are consistent with the predictions from SM processes within two standard 
deviations, and are summarized in Tables [2] and [3] Table [2] shows the number of events for 
signal, total background, and observed data, while Table [3] shows the detailed breakdown of 
backgrounds from different sources. 


Events observed in the b jet pair and Z + jet mass plane, together with the predictions from 
background processes, are plotted in Fig. [6j The b jet pair and Z + jet mass distributions are 
shown separately in Fig. [7] 
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Figure 6: Distributions of the Z + jet mass versus b jet pair mass in data (top left and right), 
and estimated background (bottom left and right). The number of signal candidate events 
is counted in the rectangular boxes defined for each signal mass hypothesis, as discussed in 
Section |3| 


Upper limits are calculated on the ratio of the measured 0° pair production cross section times 
branching fraction B ( 0° —> Zg) x £>(Z —> ££) x B(0 ] —? bb) x 2, to the theoretically expected 
cross section times branching fraction, for each search region separately. We have used as a 
benchmark the case where £>(0° —> bb) = £>(0° —> Zg) + £>(0° —> 7 g) = 0.5. The resulting 
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6 Results 



Figure 7: The b jet pair mass distributions for three different ranges of Z + jet mass in the 
electron channel (top left) and the muon channel (top right); Z + jet mass distributions for three 
different ranges of b jet pair mass in the electron channel (bottom left) and the muon channel 
(bottom right). The plotted regions correspond to three of the search regions. Predicted signal 
distributions are overlaid. The shaded band represents the statistical uncertainty combined 
with the systematic uncertainty in the simulated samples. 

limits placed on octo-triplet particle production cross section and mass can be scaled for differ¬ 
ent choices of £?(©° —> bb) and B (0° —> Zg). The results in the electron and muon channels 
are consistent and therefore combined. A modified-frequentist approach, CL S , is used to cal¬ 
culate the limits |50H5l]. For each mass point, pseudo-experiments are run for the signal plus 
background hypothesis, and for the background only hypothesis. The systematic uncertainties 
discussed in Section |5j and their correlations, are included in the limit determinations through 
a set of nuisance parameters. The results are shown, in Fig.[8j as a 95% confidence level (CL) ob¬ 
served and expected limit on a x F>(©° —> Zg) x B(Z —> ££) x £?(0° —> bb) x 2, as a function 
of 0° mass. For the case where me = 2.3m qo, this result excludes 0° masses below 623 GeV 
at 95% CL, with an expected exclusion of 639 GeV. For the case where me = 5m @ o, 0° masses 
below 426 GeV are excluded at 95% CL, with an expected exclusion of 439 GeV. These plots 
also include the theoretical predictions, with the band indicating the uncertainty in the signal 
due to the PDF uncertainty. 
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Table 2: The number of events after final selection for the signal, total background, and ob¬ 
served data, together with the statistical and systematic uncertainty For the entries where a 
single uncertainty is shown, the statistical uncertainty is negligible. 


m @ o [GeV] 

Signal 

Background 

Observed 


m G i = 2.3m 0 o 

me = 5m 0O 



Electron channel 

217 

2110±70±290 

1110±40±160 

358±15±81 

336 

304 

448±12±70 

188±6±31 

123±6±27 

115 

391 

119±3±22 

38.6±1.3±7.4 

42.3±3.0±11.0 

59 

478 

35.8±1.0±6.6 

9.45±0.32±1.90 

18.9±2.7±6.4 

24 

565 

11.7±0.3±2.4 

2.37±0.09±0.53 

8.14±1.94±3.40 

7 

652 

3.96±0.13±0.86 

0.72±0.03±0.17 

3.35±0.82±1.50 

4 

739 

1.42±0.06±0.32 

0.23±0.01±0.05 

1.02±0.31±0.51 

1 

826 

0.56±0.02±0.13 

0.08± 0.02 

0.66±0.33±0.36 

0 

913 

0.24±0.01±0.06 

0.03± 0.01 

0.31±0.14±0.19 

0 

Muon channel 

217 

2360±70±310 

1170±40±160 

348±11±74 

355 

304 

486±12±74 

185±6±31 

126±7±27 

127 

391 

118±3±20 

38.7±1.2±7.3 

44.8±3.1±11.3 

39 

478 

38.7±1.0±7.4 

9.58±0.32±1.90 

17.2±1.9±5.7 

15 

565 

13.1±0.4±2.7 

2.70±0.09±0.58 

9.52±1.40±3.80 

8 

652 

4.58±0.14±0.99 

0.74±0.03±0.17 

3.25±0.84±1.30 

5 

739 

1.63±0.07±0.36 

0.24±0.01±0.06 

1.48±0.49±0.58 

3 

826 

0.69±0.03±0.16 

0.08± 0.02 

0.32±0.17±0.14 

2 

913 

0.26±0.01±0.06 

0.03± 0.01 

0.14±0.14±0.09 

1 


19.7 fb (8 TeV) 

-.11 Trrj-rr-t^rp " 

pp -> @°@° -»(Zg)(b5), with rn,, = 2.3 m 0O 

.Expected limit 

Expected ± 1 a 
Expected ± 2 o 
—Observed limit 

si Color-Octet (with PDF uncertainty) - 



10 ° = 

= B(0° -> bb) = 0.5 
I B(0° -> Zg) + B(0 C 

t r i ii 111,11 

J 200 300 


600 700 800 900 

m 0 [GeV] 



m 0 [GeV] 


Figure 8: The 95% CL expected and observed upper limits on the cross section times branching 
fraction, crx B(©° —* Zg) x B(Z — > ££) x B(©° —> bb) x 2, as a function of 0° mass, for the 
case where me = 2.3m@o (left) and me = 5m &0 (right) with the band on the color octet the¬ 
oretical prediction indicating uncertainty in the signal yield due to PDF. These results assume 
£>(©° ->• bb) = 0.5 and B(©° Zg) + B(©° -)■ yg) = 0.5. 
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6 Results 


Table 3: The number of background events, from all sources, after final selection, together with 
the statistical and systematic uncertainty. For the entries where a single uncertainty is shown, 
the systematic uncertainty is negligible. 


m @0 [GeV] 

Z + jets 

tt and tW 

Diboson 

Electron channel 

217 

256±15±79 

91.3±2.2±18.0 

11.1±0.9±2.2 

304 

90.5±5.7±26.0 

27.5±1.2±5.3 

4.55±0.60±0.90 

391 

33.6±2.9±11.0 

6.55±0.60±1.30 

2.10±0.35±0.42 

478 

16.9±2.7±6.4 

1.24±0.25±0.24 

0.78±0.21±0.16 

565 

7.36±1.90±3.30 

0.36±0.11±0.07 

0.42±0.17±0.08 

652 

3.08±0.81±1.50 

0.04±0.03±0.01 

0.23±0.12±0.05 

739 

1.02±0.31±0.51 

<0.04 

<0.01 

826 

0.65±0.33±0.36 

<0.04 

0.01T0.01 

913 

0.27±0.13±0.19 

<0.04 

0.04±0.04±0.01 

Muon channel 

217 

244±10±72 

91.5±2.2±18.0 

12.7±1.0±2.4 

304 

92.2±7.2±27.0 

28.7±1.2±5.5 

5.57±0.59±1.10 

391 

35.1±3.0±11.0 

7.51±0.62±1.40 

2.20±0.36±0.42 

478 

15.0±1.8±5.7 

1.20±0.26±0.23 

1.00±0.24±0.19 

565 

8.89±1.40±3.80 

0.36±0.13±0.07 

0.26±0.11±0.05 

652 

3.03±0.84±1.30 

0.08±0.07±0.01 

0.14±0.05±0.03 

739 

1.27±0.48±0.58 

0.12±0.07±0.02 

0.09±0.04±0.02 

826 

0.23±0.16±0.14 

0.05±0.05±0.01 

0.05±0.03±0.01 

913 

0.14±0.14±0.09 

<0.05 

<0.05 
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7 Summary 

A search for pair production of neutral color-octet weak-triplet scalar particles (0°) has been 
performed based on processes where one 0 ° decays to a pair of b quark jets and the other to a 
Z boson plus a jet, with the Z boson decaying to a pair of electrons or muons. This analysis is 
based on data collected with the CMS experiment in proton-proton collisions at y/s = 8 TeV, cor¬ 
responding to an integrated luminosity of 19.7 fb 1 . The number of observed events is found to 
be in agreement with the standard model prediction. The CL S method is used to set a 95% confi¬ 
dence level limit on the cross section of octo-triplet particles, assuming 13(0° —> bb) = 0.5, with 
the remaining 0° branching fraction shared between Zg and yg. By comparing the theoretical 
predictions of the octo-triplet model and the observed limits, masses of 0° below 623 GeV for 
niQi = 2.3m 00 , and below 426 GeV for itiqi = 5m@o, are excluded at 95% confidence level. These 
are the first direct experimental bounds on the 0 ° production model. 
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